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We show that the optical and electronic properties of nanocrystalline silicon can be efficiently 
tuned using impurity doping. In particular, we give evidence, by means of ab-initio calculations, that 
by properly controlling the doping with either one or two atomic species, a significant modification 
of both the absorption and the emission of light can be achieved. We have considered impurities, 
either boron or phosphorous (doping) or both (codoping), located at different substitutional sites of 
silicon nanocrystals with size ranging from 1.1 nm to 1.8 nm in diameter. We have found that the 
codoped nanocrystals have the lowest impurity formation energies when the two impurities occupy 
nearest neighbor sites near the surface. In addition, such systems present band-edge states localized 
on the impurities giving rise to a red-shift of the absorption thresholds with respect to that of 
undoped nanocrystals. Our detailed theoretical analysis shows that the creation of an electron-hole 
pair due to light absorption determines a geometry distortion that in turn results in a Stokes shift 
between adsorption and emission spectra. In order to give a deeper insight in this effect, in one case 
we have calculated the absorption and emission spectra going beyond the single-particle approach 
showing the important role played by many-body effects. The entire set of results we have collected 
in this work give a strong indication that with the doping it is possible to tune the optical properties 
of silicon nanocrystals. 



PACS numbers: 73.22-f, 71.15-m, 78.55.-m; 78.20.-e. 

I. INTRODUCTION 

Bulk silicon is an indirect band-gap material emitting 
in the infrared region. Radiative lifetimes of excited car- 
rier are very long, causing a predominant de-excitation 
via fast non-radiative recombinations. Moreover, silicon 
has a significant free carrier absorption and Auger re- 
combination rates which makes the use of this material 
in optoelectronic applications very problematic. 

During the last decade, several breakthroughs have 
increased the hopes of using nanostructured silicon as 
an optical active material.^ ^ The basic idea has been 
to take advantage of the reduced dimensionality of the 
nanocrystalline phase (1-5 nm in size) where quantum 
confinement, band folding and surface effects play a cru- 
cial role.^ ^ Indeed, it has been found that Si nanocrystals 
band-gap increases with decreasing size with a lumines- 
cence external efficiency in excess of 23%.[l™ Moreover, 
optical gain has been already de monstrat ed in a large 
variety of experimental conditions) ^ * ^ * ^ * ^ * ^ Nevertheless, 
Si nanocrystals (Si-nc) still have a memory of the indi- 
rect band gap of the bulk phase and this is evidenced 
by the clearly obser ved str uctures related to momentum- 
conserving phonons-f^^^^l^This drawback can be circum- 
vented by introducing an isoelectronic impurities .fll^ In- 
deed, in a series of recent very interesting papers, Fujii 



and collaborator j-*^^ ! ^^ -^ have proved the possibility of a 
detailed control of the Si-nc photoluminescence by the 
simultaneous doping with n- and p- type impurities. In 
particular, they have showed that a (B and P) codoped 
Si-nc always has an higher photoluminescence intensity 
than that of both a single (B or P) doped and of an un- 
doped nanocrystal. Besides, under resonant excitation 
condition, the codoped samples did not exhibit struc- 
tures related to momentum-conserving phonons suggest- 
ing that, in this case, th e quasi-direct optical transitions 
are predominant li^U^lMI 

From a theoretical point of view, investigations of 
impurities in silicon nanostructures are very few when 
compared to the large number of papers reporting cal- 
culations for pure, undoped systems; moreover, most 
of them are based on semi-empirical approaches. An 
handful number of first-principles studies has been de- 
voted to quantum confinement effects in single-doped Si- 
nc.^^ 16 17 18 xhese works have basically shown that i) 
the Si-nc ionization energy is virtually size independent, 
ii) the impurity formation energy is greater for smaller 
nanocrystals and iii) impurity segregation stro ngly affect s 
the conductance properties of nanostructures.^^^^^^^^^^ 

We have recently started a systematic study of the 
electronic and optical properties of codoped Si-nc. Our 
first results^^ show that codoped Si-nc undergo a mi- 
nor structural distortion around the impurities and that 
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the formation energies are always smaller than those 
of the corresponding single-doped cases. Moreover, we 
have found that codoping reduces the band-gap with re- 
spect to the undoped nanocrystals, showing the possibil- 
ity of an i mpur ity based engineering of the Si-nc optical 
properties EHEol 

We report here a comprehensive theoretical study of 
the structural, electronic and optical properties of B and 
P simultaneously doped Si nanocrystals using ab initio 
Density Functional Theory. The paper is organized as 
follows. Section [ll] is focused on the description of the 
theoretical and computational methods, section III is 
dedicated to the discussion of our results. The results 
will be presented and discussed in the following order: i) 
structural properti es (su bsection III A), ii) formation en- 
ergies (subsection [III B ) and iii) electronic (subsection 
III C| ) and iv) optical properties (subsection |IIID| ) for 
all the considered Si nanocrystals. Concerning the op- 
tical properties we discuss both absorption and emis- 
sion spe ctra obtained within a single-particle approach 
(HID 1 ) and, in one case, with many-body methods con- 
sidering self-ene rgy corrections and the Bethe-Salpeter 
scheme (HID 2). Finally, in section IV 
our results. 



we summarize 



II. COMPUTATIONAL METHODS 

Our results are obtained within a plane-wave 
pseudopotential DFT scheme, using the Quantum- 
ESPRESSO package^^. The impurity states are 
calculated in an approximately spherical Si-nc, built by 
considering all the bulk Si atoms contained in a sphere 
(centered on a Si ion) with a diameter ranging from 
1.1 nm (SissHae) to 1.79 nm (Sii47Hioo). The surface 
dangling bonds are saturated with hydrogens. Following 
the work of Fujii et al.l^, we have located the B and 
P impurities in substitutional positions just below the 
nanocrystal surface. It is worth mentioning that this 
arrangement represents the most stable configura tion, as 
confirmed by theoretical and experimental works.^ ^^*^^ * 
Full relaxation with respect to the atomic positions has 
been allowed for both doped and undoped systems. All 
the DFT calculations are performed within the gener- 
alized gradient approximation (GGA) using Vanderbilt 
ultrasoft^^ pseudopotentials for the determination of 
both the structural and electronic properties (see Sec. 
Ill A[ Sec. [Ill B| and Sec. IIIC) whereas norm-conserving 
pseudopotential within the Local Density Approxima- 
tion (LDA) at the relaxed geometry have bee n use d 
to evaluate the optical properties (see Sec. HID). 



This choice is due to the fact that although Vanderbilt 
ultrasoft pseudopotentials allow the treatment of several 
hundreds of atoms per unit cell in the atomic relaxation 
process, the removal of the norm-conservation condition 
is a well known problem forthe calculation of the optical 
transition matrix elements^^. Each nanocrystal has been 
placed in a large supercell in order to prevent interac- 



tions between the periodic replicas (about 6 A of vacuum 
separates neighbor nanocrystals). Structural, electronic 
and optical properties, as well as the impurity formation 
energies, are investigated as a function of the size and 
for several impurity positions within the Si-nc. Starting 
from the SinH^ nanocrystal,^^ the formation energy of 
the neutral B or/and P impurities can be defined as the 
energy needed to insert one B or /and one P atom within 
the nanocrystal after removing one/two Si atoms (trans- 
ferred to the chemical reservoir, assumed to be bulk Si^ 

Ef = E (Sin-^-feBfeP/H^) — E (SinH^) 

+ (/c + /)/iSi - /c/iB - ^MP, (1) 

where E is the total energy of the system, /isi the total 
energy per atom of bulk Si and /iB(P) the total energy 
per atom of the impurity (we consider the total energy 
per atom in the tetragonal B50 crystal for B, as in Ref. 
1281 and the orthorhombic black phosphorus for P, as in 
Ref. [29]). The integers k and / can be set to either or 
1. In particular, k = 1 is the choice when a B impurity 
is present in the nanocrystal (0 otherwise) and / = 1 for 
a P impurity (0 otherwise). With this prescription, Eqjl] 
can be used for both the single doping and the codoping 
case. 

The calculations of the optical properties have been 
done both in the ground and in the first excited states, 
where the excited state corresponds to the electronic con- 
figuration in which the highest occupied single-particle 
state (HOMO) contains a hole, while the lowest unoc- 
cupied single-part icle state (LUMO) contains the corre- 
sponding elect ron.!^^^!^ It is worth pointing out that 
an undoped and relaxed Si-nc have T^^ symmetry; in the 
presence of doping, this high symmetry is generally lost 
due to the presence of the impurity atoms. Moreover, 
because of the significant differences in the charge den- 
sity of the ground and excited states, the actual atomic 
relaxations in the two cases are different. 

The nanocrystal optical response is evaluated for both 
the ground and the excited state relaxed geometries 
computing the imaginary part of the dielectric func- 
tion {€2(00)) through the Fermi golden rule. The emis- 
sion spectra have been calculated using the excited state 
atomic positions and the ground state electroni c config - 
uration (more details can be found in section III D 2 ) . 
It should be noted that although €2(00) should only be 
used for calculating the nanocrystal absorption coeffi- 
cient, it can also be used for getting a first approxima- 
tion to the emission spectra simply because the emission 
can be viewed as the time reversal of the absorption.!^ 
In other words, once the relaxed atomic positions corre- 
sponding to to an hole in the HOMO and an electron in 
the LUMO have been found, these atomic positions are 
used for the calculation of €2(00) whose main features are 
also those of the emission spectra. It is worth mention- 
ing that the photoluminescence spectra can be derived 
using the well known Van Roosbroeck and Shockle}^ 
relation which, again, involves €2(00). However, such a 



calculation requires the knowledge of the electron and 
hole populations, at the working temperature, in the 
LUMO and HOMO states respectively. The populations, 
in turn, depends on the actual dynamics in the excita- 
tion and emission processes, including the non radiative 
electron hole recombinations. In this work we have not 
considered any particular dynamics so that our emission 
spectra contains only the informations related to both 
the transition energies and the oscillator strengths. In 



TABLE I: Bond lengths (in A) around the impurity sites for 
the undoped, single doped and codoped SisrHre nanocrystal 
(diameter 1.50 nm). B and P impurities have been substitu- 
tionally located at subsurface positions (see Fig. [T]). Sis and 
Sii refer to two surface and two inner Si atoms around this 
site respectively. 
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the case of the SissBPHae codoped Si-nc, going beyond 
the single-particle approach, we have included the self- 
energy corrections by means of the GW approximation.*^ 
In a successive step, excitonic effects are include solving 
the Bethe-Salpeter equation.^ A further advantage of 
this procedure is that the inhomogeneity of the system 
is taken into account by properly including local fields 
effects. This approach, in which many-body effects are 
combined with a study of the structural distortion due to 
the impurity atoms in the excited state, allows a precise 
determination of the Stokes shift between absorption and 



emission spectra 
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III. RESULTS 



This section collects all the results we have obtained in 
the study of the structural, electronic and optical proper- 
ties of boron and phosphorus codoped silicon nanocrys- 
tals. When possible our outcomes are compared with 
available experimental results. 



A. Structural Properties 

First of all, it is interesting to look at the change in 
the nanocrystal structure induced by the presence of the 
impurities. As outlined above, the B or/and P impurity 
atoms have always been located in substitutional sites in 
the Si shell just below the surface, these positions hav- 
ing previously been shown to be the most stable ones.^ 




FIG. 1: (Color online) Relaxed structure of SissBPHre (d = 
1.50 nm). Yellow (grey) balls represent Si atoms, while the 
white (light grey) balls are the hydrogens used to saturate 
the dangling bonds. B (magenta, dark grey) and P (black) 
impurities have been located at subsurface position in substi- 
tutional sites on the nanocrystal opposite side. The relaxed 
impurity-impurity distance is Dbp = 10.60 A. 



Initially, we have considered impurities located on oppo- 
site sides of the nanocrystals, thus at the largest possible 
distance. 

Table I gives the relaxed bond lengths around the im- 
purities for SigrHye whose structure is shown in Fig{T] 
Comparing these bond lengths with those of the corre- 
sponding Si atoms in the undoped Si-nc it is clear that 
some significant relaxation occurs around the impurities. 
In all the cases the local structure has a C2v symmetry, 
with two shorter and two longer Si-impurity bonds with 
respect to the two surface and the two inner Si atoms. An 
interesting point is that the amount of relaxation around 
the impurity is directly related to the impurity valence. 
The most significant relaxation is found for the trivalent 
atom (B, 2.036 and 2.014 A with respect to 2.355 and 
2.363 A) to be compared with that of the pentavalent one 
(P, 2.294 and 2.380 A with respect to 2.355 and 2.363 A). 
Besides, it is interesting to note that in the codoped case 
the differences among the four impurity- Si bond lengths 
are always smaller than the single-doped case (the Si-B 
bonds differ of about 1.08% in the single-doped case and 
only 0.64% in the codoped case, whereas this variation 
in the case of P reduces from 3.61% to 1.54%). Thus, if 
carriers in the Si-nc are perfectly compensated by simul- 
taneous n- and p-type doping, an almost Td configuration 
is recovered. 

This tendency towards a symmetry of codoped Si- 
nc is also obtained for smaller and larger nanocrystals, 
showing that these outcomes are independent of the Si- 
nc size. Anyway, a symmetry lowering with respect to 
SigrHre is present due to the different neighborhood ex- 
perienced by the impurities. It should, in fact, be noted 
that in the case of SiasHse and Sii47Hioo the atoms in 
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TABLE II: Bond lengths (in A) around the impurity sites for 
the undoped, doped and codoped SissHse nanocrystal (d = 
1.10 nm). Substitutional B and P impurities are located at 
subsurface positions (see Fig. [2]). Sis and Sii label, respec- 
tively, the three surface and one inner Si atoms around this 
site. 



TABLE III: Bond lengths (in A) around impurity sites for 
the undoped, single and codoped Sii47Hioo nanocrystal (d = 
1.79 nm). Substitutional B and P impurities are located at 
subsurface positions (see Fig. [s]). Sis and Si^ have the same 
meaning as in Table [ill 
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FIG. 2: (Color online) Relaxed structure of the SissBPHse 
codoped nanocrystal(d =1.10 nm). Grey (yellow) balls rep- 
resent Si atoms, while the white (small grey) balls are the 
hydrogens used to saturate the dangling bonds. B (magenta, 
dark grey) and P (black) impurities have been located at sub- 
surface position in substitutional sites on opposite sides of the 
nanocrystals. The relaxed impurity-impurity distance is Dbp 
= 3.64 A 



FIG. 3: (Color online) Relaxed structure of the Sii45BPHioo 
codoped nanocrystal (diameter d = 1.79 nm). Yellow (grey) 
balls represent Si atoms, while the white (small grey) balls are 
the hydrogens used to saturate the dangling bonds. B (ma- 
genta, dark grey) and P (black) impurities have been located 
at subsurface position in substitutional sites on opposite sides 
of the nanocrystals. The relaxed impurity-impurity distance 
is Dbp = 13.59 A 
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FIG. 4: (Color online) "Phosphorus impurity path" in 
Sii45BPHioo- Atoms have the same color as in Fig. [2] The P 
atom (black) has been moved to explore several substitutional 
sites (labeled by roman numbers) from the position labeled 
by II to the positions III, V-a, V-b, Vl-a, Vl-b, Vlll-a, Vlll-b, 
IX-a, IX-b and X. The B atom (magenta, dark grey) position 
is fixed. 



repeated the calculation for the corresponding Si-Si dis- 
tance. The results are shown in Fig. [5] where, as a ref- 
erence, we also show the corresponding distances in bulk 
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IG. 6: (Color online) Formation energy for single-doped and 
codoped Si-nc. In the codoped nanocrystals the impurities are 
placed as second neighbors in the first subsurface shell (see 
text). Green triangles are related to SissHse, blue diamonds 
to SisrHre and red circles to Sii47Hioo based nanocrystals. 
The lines are a guide for eyes. 



depend on the distance between the impurities. 



B. Formation energy 

The different structural deformations occurring in the 
single-doped and codoped nanocrystals around the im- 
purity (see Tables I, II and III) have a deep influence on 
the stability of the analyzed systems. As stated in Sect. 
•j^Tp Parting from Si^H^,^ the formation energy of the 
n-^qt^al-^gr^n^ ;jj^im)u^it-^s can be defined as the en- 
ergy neeci^^"^ insert one ^ oi|/and one P atom within 
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In order 1;o clarify whiOMre the parameters that play 
an importaijit role in the aeterriination of the formation 
enetgy, we HavejMformed a sef es of total energy calcula- 
tioris considmnPi) single-dop( d and codoped nanocrys- 
tals, ii) naiRrystals of differert sizes, iii) impurities lo- 
cat^ in different sites and iv) v iriable impurity-impurity 
di^nce within a nanocrystal." 

In Fig. |6]we report the calculated formation energies 
of SissHse (diameter d= 1.10 nin), SigrHye (d= 1.50 nm) 
and Sii47Hioo (d = 1-79 nm). In the same figure, as 
a reference, we report also lih^ single-doping formation 
^Ealt^s. For the codoped case, B and P impurities have 
jd neighbors . This choice corresponds 
)o^ibley^^ta]ic(; between two subsurface 
5i33B*^6T^^e":^ig. 2| and SImsBPHioo 
I(feMtlie?f)Si^oiJft &tJ)^H^I in Fig. [4)). After a geometry 
(JSiSfii^^^Pi (^Ij^ces bet\^een B and P impurities are 
Dbp = 3.56 3.64 A and Dbp = 3.68 A for 




one /two Si atoms. 



FIG. 5: (Color online) B-P substitutional sites distances in 
Sii45BPHioo (blue circles) and the corresponding Si-Si dis- 
tances in the undoped Sii47Hioo nanocrystal (cyan squares) 
compared to the corresponding distances in bulk Si (black 
diamonds) . 



SiasBPHse, SissBPHre and SiUsBPHioo respectively. 

From Fig. [6] it is clear that the simultaneous B and 
P doping strongly reduces (of about 1 eV) the forma- 
tion energy with respect to both B and P single-doped 
cases and that this reduction is similar for Si-nc of dif- 
ferent sizes. Thus, while B or P single doping is very 
costly (in particular, the formation energy increases with 
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FIG. 7: (Color online) Formation energy of single-doped and 
codoped SisrHre and Sii47Hioo nanocrystals. Two different 
impurity-impurity distances are considered in the codoped 
nanocrystals (dashed and solid lines, larger and smaller dis- 
tance respectively, see text). Red circles refer to SisrHre, 
blue squares to Sii47Hioo- The lines are a guide for eyes. 



FIG. 8: (Color online) Formation energy as a function of the 
boron-phosphorus distance. Roman numbers label the posi- 
tions of the P atom (see Fig. |4|. The dotted and dashed 
lines connect two subsets of impurity sites in which the sur- 
rounding surface Si atoms are bounded to the same number 
of passivating H atoms. 



decreasing the nanocrystals size, in agreement with pre- 
vious calculations^^ ^^), the codoping is much easier and, 
as a good approximation, independent of the nanocrys- 
tal size. The important point here is that Si-nc can be 
more easily simultaneously doped than single-doped; this 
is due to both the charge compensation and to the minor 
structural deformation. 

It is interesting to look at the detailed dependence of 
the formation energy on the distance between the two 
impurities. In Fig. [7|we present the comparison between 
the formation energies of SigsHye and Sii45Hioo with im- 
purities placed at two different distances : 1) the previ- 
ous considered second neighbors ones and 2) the largest 
possible impurities distance (Dbp = 10.60 A and Dbp 
= 13.29 A for the SissBPHre and SImsBPHioo respec- 
tively, see Fig. [l] and Fig. [3|. We note that when the 
impurity-impurity distance is reduced, the formation en- 
ergy decreases of 0.2-0.3 eV taking negative values. This 
fact demonstrates that a stronger interaction between im- 
purities leads to a reduction in the formation energy, so 
that codoping result to be easier and more likely when 
the dopants are closer to each other. In the latter case, 
the reduction of the formation energy is almost indepen- 
dent of the nanocrystal size, as shown in Fig. [6j 

In order to investigate in more detail the dependence of 
the formation energy on the impurity-impurity distance, 
we focus our attention on the codoped Sii45BPHioo, try- 
ing to trace a "formation energy path" by progressively 
increasing the B-P distance. In this calculation we have 
kept the B atom frozen in a subsurface position while 
moving the P atom through different substitutional sites 
along the first subsurface shell, as schematically sketched 
in Fig. |4] The results of these calculations are shown in 
Fig. I 

Two interesting effects are evidenced in this figure. 
The first one is that the formation energy assumes a neg- 
ative value when the impurities are placed at distances 
smaller than 10 A, evolving towards positive values for 



larger distances. This change of sign can lead to the def- 
inition ( 
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FIG. 9: (Color online) Formation energy as a function of Si-nc 
size and impurity-impurity distance. For the three different 
Si-nc one has the following nanocrystal diameters (d) and 
impurity-impurity distances (Dbp): SissBPHse d= 1.10 nm 
and Dbp = 3.56 A (site II); SissBPHre d= 1.50 nm and Dbp 
= 2.00 A (site I), Dbp = 3.64 A (site II), Dbp = 10.60 A 
(site VIII); SImsBPHioo d= 1.79 nm and Dbp = 3.68 A (site 
II), Dbp = 4.40 A (site III), Dbp = 6.93 A (site V-a), Dbp 
= 7.91 A (site V-b), Dbp = 9.95 A (site Vl-a), Dbp = 11.32 
A (site Vl-b), Dbp = 11.30 A (site Vlll-a), Dbp = 11.93 A 
(site Vlll-b), Dbp = 12.20 A (site IX-a),DBP = 13.29 A (site 
IX-b), Dbp = 13.59 A (site X). 
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distance between impurities plays a fundamental role on 
the decrease of the formation energy. For each nanocrys- 
tal, the formation energy takes on negative values be- 
low a given distance. Moreover, the formation energy 
have a minimum value when the impurities are located 
at the minimum possible distance. Indeed, the impurity- 
impurity distance seems to play a major role with respect 
to the nanocrystals size, since the formation energy for 
similar impurity configurations are quite independent of 
the nanocrystal dimension. The small difference between 
the SigsBPHre and the SiasBPHag and Sii45BPHioo is 
due to the different neighborhood experienced by the im- 
purities in the three cases (see Tables I, II and III). 

Another relevant point is the possibility to identify two 
distinct trends for the formation energy (see Fig. 8| that 
can be related to the type of silicon cage surrounding the 
P dopant site. One can group (dotted line) together the 
cases in which the P impurity is located in the positions 
labeled II, Vl-b and X with respect to the B impurity (see 
Fig. [4|. In these positions two of the surface Si atoms 
bounded to the P impurity present two passivating H 
atoms instead of one, a situation that dominates in all 
the other configurations. A different number of capping 
H atoms influences the formation energy. 



C. Electronic Properties 

In this section we will investigate the role of codoping 
on the electronic properties of Si-nc. As in the corre- 
sponding bulk system, the insertion of impurities tends 
to modify the electronic structure. We shall show that, 
by properly controlling the doping and the size, it is pos- 
sible to modulate both the electronic structure and some 
optical features. In particular, we shall show how the 
electronic properties of the codoped nanocrystals depend 
on both the nanocrystal size and on the distance between 
the impurities. Some of the results will be discussed in 
terms of wave function localization around the impurity 
complex. 

In the single-doped cases we have already shown that 
the presence of either donor or acceptor states can con- 
siderably lower the energy gap (Eg, th e HOM O-LUMO 
energy difference) of the undoped Si-nc,'^^^ defining in 
this case the energy gap as the gap between the impu- 
rity level (partially filled considered as the HOMO) and 
the LUMO (which is empty). In these cases the par- 
tially filled HOMO level is strongly localized either on 
the B or on the P impurity. For example, in the case of 
SigeBHye the defect level is located just 0.28 eV above the 
valence band reducing the above defined energy gap from 
2.59 eV (the value for the undoped Si-nc) to 2.31 eV. In 
SigePHre the defect level is located just 0.28 eV below 
the conduction band so that the energy gap is only 0.28 
eV.^ It is interesting to note that the experimental sub- 
stitutional donor binding energy for P in bulk Si is about 
33 meV, while the experimental acceptor energy for B in 
Si is 45 meV,^ showing how, in the case of nanocrystals. 
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FIG. 10: (Color online) Calculated energy levels of the un- 
doped Sii47Hioo (left panel), codoped Sii45BPHioo with an 
impurity-impurity distance Dbp = 13.29 A (central panel), 
codoped Sii45BPHioo with an impurity-impurity distance 
Dbp = 3.68 A (right panel). The alignment has been done 
locating at the same energy the fully occupied levels with the 
same type of localization. Site X and Site II are referred to 
Fig. [4] H stands for HOMO, L for LUMO. 



the combined effects of both quantum confinement and 
weak screening tend to "transform" shallow impurities in 
"deep" centers.™!" 

The electronic properties of B- and P- codoped Si-nc 
are qualitatively and quantitatively different from those 
of either B- or P- single-doped Si-nc. Now the system is 
a semiconductor and the presence of both the impurities 
leads to a HOMO level that contains two electrons and 
to a HOMO-LUMO energy gap strongly lowered with 
respect to that of the corresponding undoped nanocrys- 
tals. Fig. [To] shows the energy levels of Sii47Hioo and 



Sii45BPHioo with the impurities located at two different 
distances. 

In one case, the impurities are placed at the largest 
possible distance (Dbp = 13.29 A) and in the other one 
at the already discussed minimum distance (Dbp = 3.68 
A) for this particular nanocrystal. From the figure it is 
evident that when impurities are at the larger distance. 
Eg is strongly reduced with respect to the corresponding 
undoped value {Eg is lowered from 2.30 eV to 1.63 eV). 
On the contrary, when the impurities are close to each 
other. Eg enlarges {Eg = 2.03 eV) although it still re- 
mains below the undoped case. We can think that when 
impurities are brought closer, the Coulomb interaction 
becomes stronger so that the energy gap becomes larger. 
Boron and phosphorus feel each other like a B-P com- 
plexes with a gap opening recalling the DFT-LDA cal- 
culated gaps of the boron phosphide bulk system: direct 
gap (F ^ F) 3.3 eV, indirect gaps (F ^ X) 2.2 eV, 
and (F A) 1.2 eV, as described in Ref. SH 

These behaviors are corroborated by the calculated 
HOMO and LUMO wave functions. Fig. [TT] shows the 
square modulus contour plots of the HOMO and LUMO 
states of the two considered Sii45BPHioo nanocrystals. 

The top panel shows the contour when the impurities 
are at a large distance while the bottom panel is that 
with the impurities at short distance. It clearly appears 
from these contours that on going from the case with well 
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FIG. 11: (Color online) From the top to the bottom: the 
HOMO (left) and LUMO (right) square modulus contour 
plots calculated for Sii45BPHioo (atom colors same as in 
Fig. |2] The impurities are located on opposite sides of the 
nanocrystal, at distance Dbp = 13.29 A (a) or as second neigh- 
bors, with an impurity-impurity distance Dbp = 3.68 A. The 
isosurfaces correspond to 10% of the maximum value. 
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oftfforlfme^f' Comparison between energy gaps of 
the undoped ^lac!c tnangles) and the codoped (red circles) 
nanotrystS^ cfe a function of the nanocrystal radius. Impuri- 
ties are located in the first shell below the surface, as second 
neighbors. The impurity-impurity distances are 3.56 A, 3.64 
A, and 3.68 A for SissBPHse, SigsBPHre, Sii45BPHioo re- 
spectively. The lines are a guide for the eye. 



increase of the impurity distance; moreover, also in this 
case we can figure out the presence of two different trends 
related to the different surface region experienced by the 
P atom in the sites II, Vl-b and X, with respect to the 
other ones (see also Fig. [sjand related discussion). Fig. 
T^pdiilt^ DUt'hiolv.' atll^a 'st' ilQ principle, it is possible to 

the_ impurity-impurity distance, 
for 3i-nc larger than those con- 
possible by codoping to obtain 
a energy gap even smaller than that of bulk Si. 
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The possibility of modula_t ng the electronic properties 
of the codoped Si-nc is also e\ ident if we keep the distance 
between the impurities constant and look at the depen- 
dence of the energy gap on the Si-nc size. Fig. [13] shows, 
for three different nanocrysfiils where the impurities are 
I^aced as second neighbors, _iow the undoped nanocrys- 
tal^sia^rgy .^^Pjis reduced in tie presence of codoping (see 
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FIG. 12: (Color online) The HOMO-LUMO energy gap as 
a function of the distance between the B and P impurities 
within the Sii45BPHioo nc. Roman numbers refer to the posi- 
tions of the P atom (see Fig. |4|. The dashed and dash-dotted 
lines connect the two subsets of impurity sites in which the 
surrounding surface Si atoms are bonded to the same number 
of passivating H atoms. 



ga^-^^f Smaller 5gSR^ry"^ls) is observed for both the 
undop^Iiiii' the codo^SSiSises. Moreover, the energy 

is " shifted towards lower ener- 
t le undoped Eq; this shift is 
sj]^ynge]| tor snjaj\er ngy^Vcr^^ Playing with both the 
nanocrystal size^and the distance between the impurities 
(Aijw interesting routes may be opened for optoelectronic 
applications. Looking at the energy gap trends in Fig. 
T2]and Fig. [13] and considering that in the codoped case 



Fujii et al.^^ found photoluminescence peaks centered in 
the 0.9-1.3 eV energy region, we may conclude that Si- 
nc playing a role in the experiment have dimensions of 
the order of few nanometers. This conclusion is consis- 
tent with the experimental outcomes^ that indicates an 
average nanocrystal diameter of about 5 nm. 



TABLE IV: HOMO-LUMO gap (Eg) for SissHse, SigrHre, 
Sii47Hioo and the corresponding codoped SissBPHse, 
SissBPHre, and Sii45BPHioo. Impurities are second neigh- 
bors. The impurity- impurity distances are 3.56 A, 3.64 A, 
and 3.68 A, for SissBPHse, SigsBPHre, and SiMsBPHioo re- 
spectively, d is the nanocrystal diameter. 



starting nc 


d (nm) 


Eg undoped (eV) 


Eg codoped (eV) 


SissHse 


1.10 


3.51 


2.86 


SigrHre 


1.50 


2.59 


2.29 


Sii47Hioo 


1.79 


2.30 


2.03 



D. Optical Properties 

The aim of this section is to investigate the mecha- 
nisms involved in the modification of the optical prop- 
erties of codoped Si nanocrystals. We present absorp- 
tion and emission spectra with a comparison between the 
IP-RPA (independent particle-random phase approxima- 
tion) spectra and the many-body ones. These last ones 
are obtained within a GW-BSE approach that takes into 
account not only the self-energy correction and the lo- 
cal field effects but also the electron-hole interaction. All 
the calculations performed are not spin-polarized. How- 
ever it should be noted that single-particle calculations 
for undoped Si-nc have been done by Franceschetti and 
Pantelides^^ within the local spin-density approximation, 
showing that the singlet-triplet splitting is significantly 
smaller than the Stokes shift. To understand the role 
of dimensionality and impurity distance and to show the 
importance of including many-body effects in the optical 
spectra, we are going to present first the result of a RPA 
independent particle optical response for various codoped 
nanocrystals different in dimensions and in impurity lo- 
cation (see Sect. 



IIIDl), and next, we will present a 
complete study of a codoped Si-nc where we go beyond 
the single-particle approach within the GW-BSE frame- 
work (see Sect. |IIID2). 



1. Absorption and emission spectra: single-particle results 

We first discuss the results related to the absorption 
spectra. Fig fM] shows a comparison between the undoped 
SigrHre and the codoped SigsBPHye IP-RPA absorption 
spectra. In this case the impurities are located at a dis- 
tance of 10.60 A. The optical response is evaluated for the 
ground state relaxed geometry computing the imaginary 
part of the dielectric function €2(00). It can be seen from 



Fig 14 that new transitions arise below the absorption on- 
set of the undoped Si-nc. In particular, we have found a 
shift of the absorption gap to lower energies with respect 
to the undoped case toghether with an enhancement of 
the intensities around 2.0 eV. These new transitions are 
due to the presence of new HOMO and LUMO state s lo- 
calized on the impurities, as described in Sec. HI C| (see 
for example Fig. pTj). The inset of Fig. 14 clarifies how 
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FIG. 14: (Color online) Comparison between the undoped 
SisrHre (solid line) and the codoped SissBPHye (dashed line) 
single-particle absorption spectra. The impurities are at a dis- 
tance of 10.60 A. Arrows indicate the energy gaps. The calcu- 
lated energy levels for the codoped nanocrystal are shown in 
the inset. A Gaussian broadening of 0.1 eV has been applied. 
H stands for HOMO, L for LUMO. 



the peak located in the 2.0-2.2 eV energy region is re- 
lated to contributions that involve the HOMO-1, HOMO 
to LUMO, LUMO + 1 transitions; it should be noted 
that for all these levels the wavefunctions are predomi- 
nantly h--'- --' Tf ^..^ ^U... 
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FIG. 15: (Color online) Single-particle absorption spectra of 
SissBPHse (black solid line), SissBPHre (red dashed line), 
and Sii45BPHioo (dash-dotted blue line). In all cases the 
impurities are second neighbors. The impurity- impurity dis- 
tances are 3.56 A, 3.64 A, and 3.68 A for SissBPHse, 
SissBPHre, Sii45BPHioo respectively. Arrows indicate the 
energy gaps. In the inset a zoomed view of the spectra onset. 
No Gaussian broadening has been applied. 
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FIG. 16: (Color online) Single-particle absorption spectra 
of Sii45BPHioo Si-nc with impurities placed at different dis- 
tances (see the inset). Arrows indicate the energy gaps. Ro- 
man numbers refer to the sites occupied by the P atom with 
respect to the B one (see Fig. [4|. No Gaussian broadening 
has been applied. 



whose diameters are 1.10 nm, 1.50 nm and 1.79 nm, re- 
spectively. In all the three nanocrystals the impurities 
are second neighbors. 

Two facts emerge from this figure. First of all, on 
increasing the nanocrystal size the absorption gap is 
strongly reduced (see arrows in Fig. 



15). Second, an 



increase of the Si-nc diameter (i.e. a decrease of the im- 
purity weight with respect to the total number of atoms) 
results in a lowering of the intensity for the transitions 
that involve the impurities. 

The role of the impurity distance on the optical re- 
sponse has been investigated following the same ap- 
proach adopted for the electronic properties in Sec. 
|III C| In Fig{T6]the single-particle absorption spectra of 
Sii45BPHioo are reported with P impurity placed on sites 
II, III, IX, and X respectively (see for comparison Fig. [4|. 

Here we observe a shift of the absorption gap to lower 
energy on increasing the distance between the impurities 
(see arrows in Fig. 16). Moreover, also the intensity is 
affected by the impurity distance. Stronger transitions 
arise when the impurities are closer whereas the intensity 
gets lower when the impurities are at a larger distances; 
the optical transitions near the band edge (indicated by 
arrows in Fig. 16) exhibit weaker oscillator strengths. 

Now we discuss the results for the emission spectra 
and for the Stokes shift between absorption and emission. 
The nanoscrystal excitation has been studied consider- 
ing the excited state as the electronic configuration in 
which the highest occupied single-particle state (HOMO) 
contains a hole (/i), while the lowest unoccupied single- 
particle state (LUMO) contains the corresponding elec- 
tron (e)^thus simulating the creation of a electron-hole 
pair.'^^'^^EiEl] Initially the system is in its ground state 
and the electronic excitation occurs with the atomic po- 
sitions fixed in this configuration. After the excitation, 
due to the change in the charge density, relaxation oc- 
curs until the atoms reach a new minimum energy due to 



Abs. (eV) 
Ems. (eV) 
A (eV) 



SissBPHs 



23 2.4 



2.77 
1.78 
0.99 
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tion and emission energy gaps (and their 



as HOMO-LUMO differences 
relaxed geometries configura- 
are obtained within the DFT 



m the^rd 
tion, fe^ecti 

singl^iijtarticlfei-gjpproach. d is Ihe nanocrystal diameter, Dbp 
is tlijS^distance between impurities and A the calculated Stokes 



SissBPHr 



2.* 



1.50 
2.00 
2.32 
2.20 
0.12 



1.50 
10.60 
1.75 
1.36 
0.39 



the presence of the electron-hole pair. The new atomic 
positions modify the electronic spectrum, implying that 
the levels involved in the emission process change. This 
model assumes that the relaxation under excitation is 
faster than the electron-hole recombination. The diflPer- 
ence between the absorption and emission energies due to 
the differen t atomic positions represents the nanocrystal 
Stokes shift.EnMMl 

The calculations have been performed for two Si-nc of 
different sizes taking, in the larger Si-nc, the impurities 
located at different distances. As shown in Table V, both 
the absorption and emission HOMO-LUMO energies are 
affected by these two parameters. With regard to the 
first parameter, we note that the Stokes shift strongly de- 
pends on the size showing a strong decrease on increasing 
the diameter of the Si-nc. This is due to the fact that 
for larger nanocrystals the excitation determines a mi- 
nor distortion of the geometry. Concerning the second 
parameter, we see that the Stokes shift tends to slightly 
increase with B-P distance although this effect is small 
if compared with the lowering due to the increase of the 
Si-nc dimensions. The comparison between these results 
and the ones previously obtained for undoped clusters 
(0.92 eV for the SiasHse-nc^^ and 0.22 eV for the SigrHre- 
ncP^ confirm that the Stokes shifts is mainly determined 
by the nanocrystals size, but that nevertheless it depend 
slightly on the presence of the impurities and also on their 
mutual distance. 

Looking at the single-particle optical spectra in Fig. [l7 
we note that the HOMO-LUMO transition in SissBPHre 
(1.75 eV, bottom panel) is almost dark when the two im- 
purities are far apart and becomes instead allowed (2.32 
eV, top panel) when their distance decreases. As dis- 
cussed before, this oscillator strength enhancement is a 
consequence of the character of the HOMO and LUMO 
states in the two cases. The emission (red-dashed lines 
in Fig. 17) spectra is red shifted with respect to the ab- 
sorptionfblack-solid lines in Fig. 17). This red shift is 
a consequence of the geometry relaxation in the excited 
state due to the excess energy necessary for promoting of 
an electron in the LUMO level. The dependence of the 
emission spectra both on the nanocrystals size (see Table 
V and Fig. 15) and on the impurities positions (see Figs. 
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FIG. 17: (Color online) Single-particle imaginary part of the 
dielectric function for the codoped SissBPHre nanocrystal in 
the ground (black-solid line) and in the excited (red-dashed 
line) geometries. B and P atoms are at the smallest possible 
distance (2.00 A, top panel) or at the largest possible distance 
(10.60 A, bottom panel) for this nanocrystal. A Gaussian 
broadening of 0.1 eV has been applied. 



Thus, the electron- jpole interac tion is here considered also 
in the emission gepi^etry. Tne heavy GW-BSE calcula- 
tion is made consiaering a large FCC supercell with a 50 
a.u. lattice parameter. The correlation part of the self- 
energy has been calculated using 10081 plane waves, 
while 49805 plane waves have been used for the exchange 
part Y^x- Then, the full excitonic Hamiltonian is diago- 
nalized considering more than 8000 transitions. 

Fig. [is] shows the calculated absorption and emis- 
sion spectra fully including the many-body effects. The 
electron-hole interaction yields significant variations with 
respect to the single-particle spectra (see for a compari- 
son Fig. 



17), with an important transfer of the oscillator 



16 and 17) reveals once more the possibility of tuning the 



optical response of silicon nanocrystals. 



2. Absorption and Emission Spectra: Many-body effects 

In order to give a complete description, within the 
many-body framework, of the codoped Si-nc response to 
an optical excitation, we consider both the self-energy 
corrections by means of the GW method^^ to obtain the 
quasiparticle energies and the excitonic effects through 
the solution of the Bethe-Salpeter equation. The effect 
of local fields is also included, to take into account the 
inhomogeneity of the systems. 

Since the GW-BSE calculation^^ are very computing 
demanding, we have only considered the smaller codoped 
nanocrystal SiaaBPHae (see Fig. [2|. In this particular 
cluster, we found that Local fields effects are, although 
not negligible, of minor importance with respect to GW 
and excitonic effects. It is anyway essential to include 
all of them (LF and many-body) in order to get the final 
converged spectrum shown in Fig. [Tsj 

In order to carry out emission spectra calculations, we 
use the excited state geometry and the ground state elec- 
tronic configuration. As already noted before, in this 
case €2(00) corresponds to an absorption spectrum in a 



strength to the low energy side. Moreover, in the emis- 
sion spectrum the rich structure of states characterized, 
in the low energy side, by the presence of excitons with 
largely different oscillator strengths, determines excitonic 
gaps well below the optical absorption onset. Thus the 
calculated emission spectrum results to be red shifted to 
lower energy with respect to the absorption one. This 
energy difference between emission and absorption, the 
Stokes shift, can be lead back to the relaxation of the 
Si-nc after the excitation process. 

The new important features that appear in the 
emission many-body spectra are related to the pres- 
ence of both B and P impurities as showed by Fig. 
19, which gives the real-space probability distribution 
|'0eccc(^e5 for the bound exciton as a function of the 
electron position when the hole is fixed in a given 
position. In this case the hole is fixed on the boron atom 
and we see that the bound exciton is mainly localized 
around the phosphorus atom. 

From Table VI, it can be seen that the single-particle 
DFT results strongly underestimate the absorption and 
emission edge with respect to the GW+BSE calculation, 
in which the excitonic effect are taken exactly into ac- 
count. This means that, in this case, the cancellation be- 
tween GW gap opening (which gives the electronic gap) 
and BSE gap shrinking (which originates the excitonic 
gap) is only partial. 
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eV is almost dark while an impor- 
i^^kTs evident at about 2.75 eV (see Fig. 



GfTspectrum , 
tant excitoni i 

agaij^tif^c^Bfted with respect to the first absorption 
peak.-^Ss expectedj^hat comes out is the importance 
fully taking ing^ccount the many-body aspect of the 
PSlem ii^^i^e^^^^tercome the limits of the single- 
particle ^(pproach. 

IV. CONCLUSIONS 



FIG. 19: (Color online) Excitonic wave function of SissBPHse 
(atom colors as in Fig. [2]). The grey isosurface represents the 
probability distribution^exc(re, r^)^ 
hole fixed on the B impurity. 



of the electron with the 



The difference between the GW electronic gap and the 
GW+BSE optical excitonic gap gives the exciton binding 
energy E5. We note the presence of exciton binding en- 
ergies as big as 2.2 eV, which are very large if compared 
with bulk Si (~ 15 meV) or with carbon nanotube^^^^ 
where E5 ~]^eV, but similar to those calculated for un- 
doped Si-n(j^ of similar size and for Si and Ge small 
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nanowires.i 



TABLE VL Absorption and Emission energies calculated 
as HOMO-LUMO energy difference within the single-particle 
DFT, the many-body GW and the GW+BSE approaches. A 
is the calculated Stokes shift between absorption and emission 
energy gap. The 2.20 eV energy corresponds to an almost 
dark transition. 



SissBPHse 


DFT 


GW 


GW+BSE 


Abs. (eV) 


2.80 


5.52 


3.35 


Ems. (eV) 


1.79 


4.37 


2.20 


A (eV) 


1.01 


1.15 


1.15 



The differences between full many-body calculations 
and single-particle results are of 0.55 eV and 0.41 eV for 
absorption and emission energy gaps respectively, and of 
0.14 eV between the two Stokes shifts. It is interesting 
to note that the HOMO-LUMO transition in the emis- 



he structural, electronic and optical properties of 
T rflbocrystals codoped with B and P impurities have 
been s'^died also going beyond the single-particle ap- 
proach. We have considered Si-nc of different size and 
with the impurities located at different distances. We 
show that codoping is always energetically favored with 
respect to simple B- or P-doping and that the two im- 
purities tend to occupy nearest neighbor sites near the 
surface rather than other positions inside the nanocrys- 
tal itself. Our results demonstrate that the codoped 
nanocrystals present valence and conduction band-edge 
states which are localized on the two impurities respec- 
tively and band-gaps always lower in energy with respect 
to that of undoped Si nanocrystals. Besides, the elec- 
tronic properties show a dependence on both nanocrys- 
tal size and impurity-impurity distance. The impurity 
located band-edge states originate absorption thresholds 
in the visible region which are shifted lower in energy 
with respect to the undoped case. Moreover, the emission 
spectra show a Stokes shift with respect to the absorp- 
tion which is due to the structural relaxation after the 
creation of the electron-hole pair. Our results make evi- 
dent the presence of electronic quasi-direct optical transi- 
tions between donor and acceptor states that can help to 
understand the experimental outcomes and makes it pos- 
sible to engineer the absorption and emission properties 
of Si nanocrystals. 
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